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ABSTRACT 

Supernova remnant (SNR) Kes 69 is morphologically characterized by brightened radio, infrared, 
and X-ray emission on the southeastern rim, with the 1720 MHz OH masers detected in the north- 
eastern and southeastern regions at various local standard rest (LSR) velocities. We have performed a 
millimeter observation in CO and HCO + lines toward Kes 69. From the northeastern compact maser 
region, 12 CO and 13 CO emission's peaks around 65 kms -1 and 85 kms -1 , which are consistent with 
the masers' LSR velocities, are detected. In the southeast, a molecular ( 12 CO) arc is revealed at 77- 
86 kms -1 , well coincident with the partial SNR shell detected in the radio continuum and mid-infrared 
observations. An 85 kms -1 HCO + emission is found to arise from a radio peak on the shell. Both 
the molecular arc and the HCO + emission at ~ 85 kms -1 seem to be consistent with the presence 
of extended OH masers along the southeastern boundary of Kes 69. The morphology correspondence 
between the CO arc and other band emission of the Kes 69 shell provides strong evidence for the asso- 
ciation between SNR Kes 69 and the ~ 85 kms -1 component of molecular gas. The multiwavelength 
emissions along the southeastern shell can be accounted for by the impact of the SNR shock on a 
dense, clumpy patch of molecular gas. This pre-existing gas is likely to be a part of the cooled debris 
of the material swept up by the progenitor's stellar wind. The association of SNR Kes 69 with the 
molecular cloud at the systemic velocity of ~ 85 kms -1 enables us to place the SNR at a kinematic 
distance of 5.2 kpc. 

Subject headings: ISM: individual (Kes 69, G21.8— 0.6) - ISM: molecules - supernova remnants 



1. INTRODUCTION 

The progenitors of core-collapse supernovae are most 
probably formed in giant molecular clouds (MCs). Due 
to the short lifetime, they are not far away from their 
matrices when they explode. Therefore it is common 
that the supernova remnants (SNRs) are located in the 
vicinity of MCs and may encounter them in evolution. 
The association with MCs often results in irregular mor- 
phology of the SNRs in multiwavelengths, which indi- 
cates sophisticated shock interaction with the inhomo- 
geneous environmental medium. About 20 SNRs have 
been discovered to be physically interacting with ambi- 
ent molecular gas based on the detection of the 1720 MHz 
OH masers (Frail et al. 1996), which are believed to be 
a tracer of the shock interaction with MCs (Lockett et 
al. 1999; Frail & Mitchell 1998; Wardle & Yusef-Zadeh 
2002). 

SNR Kesteven 69 is thought to be probably associ- 
ated with MCs, because of the OH masers detected to- 
ward this remnant; however, the masers are found at 
various local standard rest (LSR) velocities and at var- 
ious projected locations. Green et al. (1997) detected a 
compact OH maser using the Very Large Array (VLA) 
and the Australia Telescope Compact Array (ATCA) at 
Vlsr = 69.3 kms -1 , which is located projectionally in 
the northeastern part of the remnant. Recently Hewitt, 
Yusef-Zadeh, & Wardle (2008) not only found that this 
compact maser also has faint emission at 85 kms -1 , but 
also detected extended OH maser emission at the veloc- 
ity of ~ 85 kms -1 with the Green Bank Telescope ob- 
servation and the VLA archival observation toward the 
southern bright radio shell. The different LSR systemic 



velocities imply the MCs at different distances that may 
be impacted by the SNR shock wave. Therefore inves- 
tigation is needed to clarify at which systemic velocity 
the maser emission is the product of the Kes 69 shock 
interaction. 

Kes 69 has an irregular X-ray morphology, as observed 
by ROSAT and Einstein, inside an incomplete radio shell 
(Seward 1990; Yusef-Zadeh et al. 2003). The Spitzer In- 
frared Array Camera (IRAC) mid-infrared observation 
toward Kes 69 shows an arc at 4.5/xm in the same loca- 
tion as the southeastern radio shell (Reach et al. 2006). 
The extended OH emission along the southern incom- 
plete radio and Infrared (IR) shell hints an interaction 
of the SNR with the dense molecular gas in the south. 
If the SNR/MC association is established, a big progress 
can be made toward resolving the open question of the 
disparate velocity maser components and the distance to 
Kes 69 can also be determined. 

Motivated by the supposed association of Kes 69 with 
molecular gas, we have performed millimeter CO and 
HCO + observations toward this remnant. The observa- 
tions and results are described in §2 and §3, and the 
conclusion is summarized in §4. 

2. OBSERVATION AND DATA REDUCTION 

The observations of millimeter molecular emissions to- 
ward SNR Kes 69 were made in two epoches during 2006 
November-2007 January and 2007 October-November 
with the 13.7 m millimeter-wavelength telescope of the 
Purple Mountain Observatory at Delingha (hereafter 
PMOD). An SIS receiver was used to simultaneously ob- 
serve the 12 CO (J=l-0), 13 CO (J=l-0), and C 18 (J=l 
0) lines. We mapped a 25' x 25' area that contains the 
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full extent of SNR Kes 69 via raster-scan mapping with 
the grid spacing of 30"-60". The main-beam efficiency 
in the observing epoch was 67% and elevation calibration 



was performed. The typical system temperature was 
around 140-280 K. Three Acousto-Optical spectrome- 
ters (AOS) were used as the back end, and the cor- 
responding spectral coverages were 145 MHz for 12 CO 
(J=l-0) and 43 MHz for both 13 CO (J=l-0) and C 18 
(J=l— 0), all divided into 1024 channels. The observed 
velocity ranges were — 120-(+260) kms -1 for 12 CO, 11- 
127kms -1 for 13 CO (J=l-0), and 10-128kms -1 for 
C 18 (J=l-0). We also chose two points for the HCO+ 
(J=l— 0) line observation with long-time integration [122 
minutes for point (18 h 32 m 50 s , -10°12'42") and 134 min- 
utes for (18 h 33 m 10 s , -10°12'42")], with the velocity 
range —3 to +141 kms -1 . The baseline subtraction was 
performed with low-order polynomial fit. In Table [1] we 
list some observational parameters, such as the frequen- 
cies, the average rms noises of all final spectra, the chan- 
nel separations, and the FWHM. All data were reduced 
using the GILDAS/CLASS package 1 . 

The Hi line emission data of the archival VLA Galac- 
tic Plane Survey (VGPS; Stil et al. 2006) were also pro- 
cessed. The observations were made by the VLA and 
the Green Bank Telescope, presented with the beam-size 
of 1', the velocity resolution of 1.56 kins" 1 , and the rms 
noise of 2 K per 0.824 kms -1 channel. The processed 
Spitzer IRAC (Fazio et al. 2004) Basic Calibrated Data 
were used here (PID: 146, PI: Ed Churchwell), which 
are available in the archive of the Spitzer Science Center 
and include flat-field correction, linearity and flux cal- 
ibrations, and dark subtraction. The final mosaic 4.5 
fim IR image was produced by further processing with 
the custom IDL software (Huang et al. 2004). We also 
used the archived Spitzer Multiband Imaging Photome- 
ter (MIPS; Rieke et al. 2004) Post Basic Calibrated Data 
to present the mid-IR 24 /im image, which were obtained 
from the Micron Survey of the Inner Galactic Disk Pro- 
gram (PID: 20597, PI: Sean Carey). The 1.4 GHz radio 
continuum emission data were obtained from the NRAO 
VLA Sky Survey (NVSS; Condon et al. 1998). 

3. RESULTS 

3.1. The CO and HCO + Emissions 

Figure Q] shows three CO spectra toward the center 
(18 h 32 m 50 s .0, -10°06'42") of the mapping area. There 
are several velocity components in the velocity range 0- 
120 kms -1 in the 12 CO (J=l-0) spectrum, and no 12 CO 
emission is detected out of this range in the whole map- 
ping area. The 12 CO emission peaks appear in the inter- 
vals - 2-12, 40-60, 60-76, 76-100, and 107-120 kms" 1 . 
We have found no morphological correspondence between 
the CO emission in the velocity intervals 2-12 kms -1 
and 40-60 kms -1 and the emission from Kes 69 in any 
other waveband. The molecular component at around 
110 kms -1 also shows no morphological correspondence, 
which is near the tangent point in this direction. Apply- 
ing the rotation curve of Clemens (1985) together with 
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Fig. 1. — The spectra of the central point 

(18 h 32 m 50 s .O, -10°06'42") of the field of view in our obser- 
vation toward Kes 69: dashed line for 12 CO, thick solid line for 
13 CO, and thin solid line for C ls O. The S/N is highest here. The 
small dip near 33 kms" 1 results from the weak emission at the 
reference position, which is far from the concerned velocities. 

R = 8.0 kpc (Reid 1993) and V = 220 km s -1 , the tan- 
gent point in this direction is at 7.4 kpc at 113.7kms -1 . 
For the 12 CO emission, we focused our analysis on the 
velocity intervals 60-76 kms -1 and 76-100 kms -1 , which 
respectively cover 69.3 kms -1 and 85 kms -1 at which the 
OH masers were detected. The 13 CO (J=l-0) emission is 
prominent at around 54 kms -1 and 82 kms -1 , at which 
the 12 CO emission is strong. The intensity of the C 18 
(J=l-0) emission is too weak to examine its spatial dis- 
tribution. 

We have produced 12 CO intensity maps with inter- 
val lkms -1 between 60 and 76 kms -1 (Fig. [2|) to ex- 
amine the molecular gas around 69.3 kms -1 , at which 
the compact OH mascr is seen. In the northeastern 
compact maser region, some faint diffuse 12 CO emis- 
sion are present at 60-70 kms -1 . A distinct cloud at 
69-76 kms -1 is seen in the west. 

We have also produced 12 CO intensity maps over the 
velocity range of 76-92 kms -1 with interval lkms -1 
(Fig. [3]). By comparison with the radio continuum emis- 
sion, an arc of the molecular gas at 77-86 kms -1 is 
strikingly coincident along the southeastern rim. This 
arc is clearly seen in the close-up intensity images at 
80-81 kms -1 (Fig. gl) and 79-82 kms -1 (Fig. [5]). In the 
northwest (see Figs. |4] and [5|), there is another section 
of molecular arc, which can even be discerned in the 
wider range 79-86 kms -1 (Fig[3]). The two sections of 
molecular arcs are seen in similar velocity range and can 
be threaded with a circle of angular radius 8'. 7. Some 
bright, complicated 12 CO features are present in the 
field of view at 80-88 kms -1 . There may be a contri- 
bution from the Hll region G21. 902-0. 368 in the north- 
west, where the recombination line at 79.5 kms -1 has 
been detected (Lockman 1989). There is also a molecular 
cloudlet at 87-91 kms -1 , coincident with the strongest 
southern radio peak, at the western end of the incom- 
plete radio shell. If this small cloud is associated with 
the SNR, the radio peak can be accounted for by the im- 
pact of the remnant shock on it, because of the drastic 
shock deceleration and hence the magnetic field compres- 
sion and amplification. 



1 http: / /www.iram.fr /IRAMFR/GILDAS 
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TABLE 1 
Observational parameters 



Line 


frequency (GHz) 


noise a (K) 


Av h (km s -1 ) 


FWHM (") 


12CO (J=l-0) 


115.271204 


0.29 


0.37 


60 


13 CO (J=l-0) 


110.201353 


0.24 


0.11 


60 


C ls O (J=l-0) 


109.782183 


0.20 


0.12 


60 


HCO+ (J=l-0) 


89.188526 


0.04 


0.14 


78 



a Average rms noise of all final spectra; b Channel separation. 
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Fig. 2. — 12 CO J=l— emission maps integrated each 1 km s — 1 (smoothed to an angular resolution of 0'.4 by linear interpolation), with 
the NVSS 1.4 GHz radio continuum emission in contours (contour levels are 20, 49, and 107 mjy beam -1 ). Central velocities are indicated 
in each image. The rms noise of each map is ~ 0.18 K km s" 1 and the beam size is 1'. The diamond denotes the location of the compact 
maser at 69.5kms — 1 . 



The CO spectra from the northeastern compact OH 
maser region and the southeastern shell region are shown 
in Figure [S] (the two regions for CO spectrum extraction 
are shown in Figure [4j. For the former region, the 12 CO 
and 13 CO emissions peak at both ~ 65 and 85 kms -1 
(both with a signal-to-noise ratio (S/N)>3), essentially 
in agreement with the velocities at which the single maser 
arises (Green et al. 1997; Hewitt et al. 2008). The 12 CO 
and 13 CO lines at ~ 65 kms -1 have complicated pro- 
files, each of which can be phenomenologically fitted with 
a combination of at least three Gaussians. The lines 
at 85 kms -1 are slightly broadened in the blue wings. 
For the latter region, there are 12 CO and 13 CO lines at 
85 km s" 1 , with blue wings broadened (although the wing 
may be contaminated by the 72 kms -1 peaks). The pro- 
file of each line is divided into a Gaussian at 85 kms -1 
and a residual part in the blue wing, with the fitted and 
derived parameters summarized in Table [2j 



In the derivation, we used two methods to estimate 
the H2 column density and molecular mass. In the 
first method, the H2 column density is estimated by the 
use of the conversion factor A(H 2 )/H / ( 12 CO)« 1.8 x 
10 20 cm -2 K -1 kirT^s (Dame et al. 2001). In the sec- 
ond method, we assume local thermodynamical equi- 
librium (LTE) for the gas and optically thick condi- 
tion for the 12 CO (J=l-0) line and use the relation 
A(H 2 ) w 7 x 10 5 A( 13 CO) (Frerking et al. 1982). 

There seems to be a weak feature at 82 kms -1 in the 
13 CO line profile of the southeastern region, which could 
be either a broadened part of the 85 kms -1 component or 
a chance coincidence of an irrelevant, unperturbed cloud 
in the line of sight. We performed a multiple Gaussian 
fit incorporating this feature and derived a molecular gas 
mass for this feature, which is about two orders lower 
than the virial mass, indicating that the cloud is per- 
turbed. In view of this and the clear shell-like appear- 
ance at 80-82 kms -1 , we treat it simply as a small part 
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Fig. 3. — The same as for Figure [2] but different velocity. 



TABLE 2 

Fitted and derived parameters for the MCs around 85kms _1 in the southeastern region 3 



Gaussian Components 


Line 


Center (km s I ) 


FWHM (km s" 1 ) 


Tp ca k (K) 


W (K km s- 1 ) 


12 C (J=l-0) 
13 CO (J=l-0) 


84.7 
84.9 


5.2 
2.4 


2.4 
0.5 


13.4 
1.2 


Molecular Gas Parameters 




N(H 2 ) (10 21 cm~ 2 ) b 


Af(M ) b 


T cx (K) c 


r( 13 CO) d 


Gaussian components: 
Residual part' : 
The wholes . 


2.4 / 1.6 

2.2 

4.6 


5.8 x 10 3 d 2 2 / 3.8 x 10 3 d 2 2 ° 

5.2xl0 3 d| 2 e 

1.1 x 10 4 d| 2 c 


7.8 


0.2 



a The region is denned in Figure|4] b See text for the two estimating methods. The excitation temperature calculated from the maximum 
12 CO(J=1-0) emission point in the region. d The optical depth of the 13 CO (J=l-0) line. c d 5 . 2 = ci/(5.2kpc) (see fl3~3t . f Determined 
by subtracting the Gaussian components centered at 72kms _1 and 85kms _1 after the multiple Gaussian fitting for the 12 CO (J=l-0) 
emission in the velocity range 70— o^kms" 1 ^ Combination of the Gaussian components around 85kms _1 and the residual part. 



in the broadened wing of the 85kms _1 gas. 

In the long-time pointing observation toward the 
southern strongest radio peak or the cloudlet mentioned 
above, which targets at (18 h 32 m 50 s .O, -10°12'42"), no 
HCO + signal was detected. However, HCO + emission 
was detected in the pointing observation toward another 
radio peak at (18 h 33 m 10 s .O, -10°12'42") on the south- 
eastern shell. The 12 CO, 13 CO, and HCO+ spectra of 
this point are shown in Figure [71 The HCO + emission 
appears to be prominent only at ~ 85kms _1 , at which 
both the 12 CO and 13 CO emission peak, too. We cal- 
culated the column density of HCO + assuming that the 
HCO + emission is in LTE with the same excitation tem- 
perature as the 12 CO emission (~ 7.8 K) and is opti- 
cally thin. For comparison with the CO column density, 



we smoothed the CO (J=l-0) data to the same angu- 
lar resolution as the HCO + observation, and the results 
are listed in Table [3] Considering the 12 CO/ 13 CO ratio 
(30-70) of the general interstellar medium (Langer 1992 
and references therein), the HCO + /CO abundance ratio 
is ~ 5xl0~ 6 . It is somewhat less than that of undis- 
turbed cold MCs such as TMC-1 or L134N (1 x 10~ 4 ; 
Ohishi, Irvine, & Kaifu 1992). The HCO+ abundance is 
expected to be reduced in a slow, nondissociating shock, 
unless some enhanced ionization is present (Iglesias & 
Silk 1978; Elitzur 1983). We note that the extended 
maser emission on the southern rim is found at the same 
LSR velocity (Hewitt et al. 2008) , which implies that the 
dense molecular gas at this velocity is disturbed by C- 
type shock. For the origin of HCO + emission, though, 



DISCOVERY OF MOLECULAR SHELLS ASSOCIATED WITH SNRs. I. KES 69 



TABLE 3 

Fitted and derived parameters for the HCO+ emitting point (18 h 33 m 10 s .O, — 10°12'42") 



Gaussian Components 


Line 


Center (km s i ) 


FWHM (km s -i ) 


Tpeak (K) 


W (K km s -i ) 


12 CO (J=l-0) 


85.0 


3.4 


3.9 


14.3 


13CO (J=l-0) 


84.9 


2.4 


0.5 


1.2 


HCO+ (J=l-0) 


85.0 


2.2 


0.2 


0.5 


Molecular Gas Parameters 


AT( 1J CO) (10 i!) cm- 2 ) 


7V(HCO+) (10 11 enr 


- 2 ) r( ia CO) 






1.3 


3.5 


0.13 
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Fig. 4. — 12 CO (J=l— 0) intensity map in the velocity interval 
80—81 kms -1 . The intensity contours (smoothed to a resolution of 
0'.24 by interpolation) are at linear scale levels from 1.67 to 5.57 
Kkms" 1 in steps of 0.6 Kkms" 1 , and the thick contour denotes 
the half-maximum intensity (2.79 Kkms -1 ). The rms noise is 
~ 0.18 Kkms -1 . The dashed circle is plotted roughly running 
through the molecular arcs in the southeast and the northwest, 
with an angular radius of 8'. 7. The diamond denotes the location 
of the compact OH maser at 69.3 km s -1 . Two regions for CO 
spectrum extraction that overlap the northeastern OH maser and 
southeastern CO arc are also shown. 

further observation is needed to make a decision. 

3.2. Morphological correlation in multiwavelengths 

Combination of the multiwavelength observations to- 
ward Kes 69 (Fig. [5]) shows some interesting morpholog- 
ical correlation. 

First, the southeastern 77-86 km s" 1 molecular shell is 
not only coincident well with the 1.4 GHz radio shell (as 
pointed out above) , but also (as seen in Fig. [8|) with the 
mid-IR "ridge" described in Reach et al. (2006). This 
mid-IR ridge is prominent in 4.5/^m, with a relatively 
high brightness in 5.8/im. A part of the 24 /an floclike 
emission in the southeast is also aligned with the molec- 
ular shell (Fig. [5]). The presence of the molecular shell is 
consistent with the HCO+ and extended 1720 MHz OH 
maser emission at 85 km s -1 along the southern rim. This 
correspondence in morphological features is most prob- 
ably another signature of the SNR-MC interaction, in 
addition to the 1720 MHz OH masers detected in Kes 69. 

Second, as described above, both the sections of molec- 
ular arcs at ~ 77-86 km s" 1 in the southeast and the 
northwest seem to be aligned along a circle. A faint 




Fig. 5. — The 12 CO (J=l-0) intensity map in the velocity in- 
terval 79— 82kms -1 (above 5<r) coded in green, overlaid with the 
ROSAT PSPC X-ray contours in white, and the NVSS radio con- 
tinuum emission in red (the contour levels are 10, 41, and 172 mjy 
beam -1 ). The diamond denotes the location of the compact OH 
maser at 69.3kms — . The dashed circle is plotted roughly running 
through the molecular arcs in the southeast and the northwest, 
with an angular radius 8'. 7. 



1.4 GHz bar in the north (highlighted by a red contour 
in Fig. [5]) is also along this circle. Moreover, as noted 
by Yuscf-Zadeh et al. (2003), there is a faint radio shell 
in the northwest in the low-resolution 330 MHz image of 
Kes 69 (Kassim 1992), which is confused with the Hll re- 
gion G21.902-0.368. We note that the peak of this patch 
of radio emission roughly coincides with the northwest- 
ern section of the molecular arc. 

These morphological correlations demonstrate that 
Kes 69 is associated with the giant MC at the systemic 
velocity of ~ 85 kms -1 . This association is strengthened 
by the detection of the extended and compact 1720 MHz 
OH maser emission and the HCO + emission at 85 kms -1 
from the SNR region. 

The southeastern half of the circle is roughly X-ray 
bright, while the other half appears to be X-ray faint but 
covered by the 79-82 kms" 1 molecular gas. This bright- 
ness anticorrelation between the X-ray and CO emission, 
however, does not imply that the molecular gas in the 
northwestern half obscures the X-rays. We estimate the 
hydrogen column density of this gas (in the interval 79- 
82kms~ 1 ) to be A H ~ 1 x 10 21 cm -3 . If there were 
an X-ray emitting gas, with similar properties to that 
observed in the southeast (temperature kT x ~ 1.6 keV 
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Fig. 6. — The CO spectra from the two regions that are defined 
in Figure [4] The upper panel is for the northeastern region cov- 
ering the compact maser and the lower panel for the southeastern 
shell region. The dashed lines stand for the 12 CO emission, the 
thick solid lines for 13 CO, and the thin solid lines for C ls O. The 
sinusoidal pattern above the noise limit that appears in the 12 CO 
spectra for negative LSR velocities originates from the nonlinear 
response of the AOS system. 




Velocity (km/s) 



Fig. 7.— The spectra of the point (18 h 33 m 10 s .O, -10°12'42") 
on the southeastern rim of Kes 69: the dashed line for 12 CO, the 
thin solid line for two times the intensity of 13 CO, and the thick 
solid line for five times the intensity of HCO+ . The noise level for 
the HCO + emission is given in Table [l] together with those for the 
CO lines. 

and intervening hydrogen column 2.4 x 10 22 cm" 2 ; Yusef- 
Zadeh et al. 2003), behind this northwestern molecular 
gas, then an extra extinction by the Nn ~ 1 x 10 21 cm~ 3 
gas would only cause an ~ 8% decrease in the 0.5-2 keV 




Fig. 8. — Spitzer 4.5 fim (upper panel) and 24 (im (bottom panel) 
grey scale images overlaid with the 80-81 km s — 1 12 CO (J=l— 0) 
intensity contours (the same as in Fig. [J}. 

X-ray flux and this X-ray emission could have been ob- 
servable. Thus we conclude that the X-ray faintness in 
the northwest is not caused by the absorption of the dif- 
fuse 79-82 kms" 1 molecular gas, whether this gas is con- 
nected to the ~ 85 kms -1 cloud or not. 

3.3. The Dynamics 

The association between SNR Kes 69 and the MC at 
systemic velocity 85 kms -1 facilitates a convincing deter- 
mination of the kinematic distance to the remnant. The 
systemic velocity 85 kms -1 is suggestive of two candidate 
kinematic distances to the SNR/MC association, 5.2 kpc 
and 9.6 kpc. The choice can be made with the aid of 
the Hi absorption along the line of sight. Following the 
method used in Tian, Leahy, & Wang (2007) and Tian & 
Leahy (2008a), we produced three Hi spectra of regions 
along the southeastern shell of Kes 69, as shown in Fig- 
ure [9] Distinct absorption features appear at 4, 6, 15, 
40, 48, 67, 70, 82, and 88 kms" 1 . The 12 CO components 
along the shell around 5 and 50 kms -1 (Fig. [6]) can, due 
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Fig. 9. — Left-top panel: HI emission map of Kes 69 from a single channel at 74.39 km s -1 , overlaid with NVSS radio continuum 
contours (with levels at 20, 49, 78, 107, 136 and 165 mjy beam ). The three pairs of regions are defined to extract the HI spectra. Other 
three panels: the source and background HI emission spectra and the absorption spectra extracted from the defined regions 1, 2, and 3, 
respectively. 



to the corresponding absorption features in the Hi spec- 
tra, can be related to the chance-coincident foreground 
gas. No Hi absorption features are present around the 
tangent point velocity 113.7 kms^ 1 , which indicates that 
the SNR is in front of the tangent point (at 7.4 kpc). 
Hence, the distance to the SNR is d = 5.2 kpc. We note 
that a similar estimate has in the meantime been given 
by Tian & Leahy (2008b). 

The southeastern molecular arc is revealed above to be 
coincident with the radio and IR shell and move at a ve- 
locity v m of order lOkms -1 (matching the blucshifted 
line broadening which reflects the velocity component 
in the line of sight). The broadened blue wing of the 
~ 85 kms -1 line profile of the molecular arc implies that 
the giant molecular cloud may have suffered a perturba- 
tion from the rear side. Three scenarios regarding the 
dynamical relation between the molecular arc and the 
SNR arc discussed below. 

First, the arc is likely to be a flake of molecular gas 



that is shocked by the slow transmitted cloud shock af- 
ter the SNR blast wave hits the molecular cloud. In 
this case, the radio continuum emission may arise from 
the SNR shock that is blocked by dense cloud (Frail 
& Mitchell 1998) or the blast wave that propagates in 
the intercloud medium (Blandford & Cowie 1982), and 
the thermal X-ray emission may be ascribed to the hot 
gas behind the shocked molecular gas or just behind the 
blast wave. Thus there can be a crude pressure bal- 
ance between the cloud shock and the X-ray emitting hot 
gas (Zel'dovich & Raizer 1967; McKee & Cowie 1975): 
nov% ~ rimV^, where n m denotes the number density 
of the hydrogen atoms ahead of the cloud shock, no the 
density of the undisturbed intercloud medium, and Vb the 
velocity of the blast wave. Velocity Vb is related to the 
postshock temperature kT x as Vb = [16kT x / (3/Itoh)] 1 / 2 , 
where tor is the hydrogen atom mass and \x = 0.61 is 
the average atomic weight. Adopting kT x ~ 1.6 keV 
from the ROSAT X-ray observation (Yuscf-Zadch ct 
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al. 2003), we have v\, ~ 1.2 x 10 3 kms _1 and n m ~ 
1.3 x 10 3 (no/0.1 cm _3 )(w m /10kms _1 ) -2 cm -3 , implying 
that the blast wave hits a very dense matter in the molec- 
ular arc (here no is assumed to be similar to the mean 
density of the X-ray emitting gas, ~ 0.1cm -3 , as ob- 
tained from a reproduced ROSAT X-ray spectral anal- 
ysis 2 ). In this scenario, assuming an adiabatic expan- 
sion and adopting the SNR extent represented by the 
circle of radius ~ 8'. 7 (gSHJ) or r s ~ 13d 5 . 2 pc, the SNR's 
age is estimated as t — (2r s )/(5«(,) ~ 4.4e?5.2 kyr and 
the explosion energy is E = (25/4£)(1.4rionJH)vf rf ~ 
6 x 10 50 (no/0.1cm -3 )df 2 erg (where £ = 2.026). 

The 4.5/im and 5.8/^m mid-IR emissions are suggested 
to be likely dominated by lines of shocked gas (Reach et 
al. 2006). They note that the 5.8/im emission is relatively 
strong, but it is not clear whether H2 or [Fe II] lines (as 
likely mechanisms) are responsible for this emission, un- 
like the case of IC 443, in which the [Fe II] and H2 lines 
are very clearly segregated to the northern and south- 
ern regions of the remnant, respectively (see Rho et al. 
2001). The H2 emission seems to be consistent with the 
slow molecular shock, while [Fe II] emission could not 
be ruled out with the present observations if the SNR 
shocks propagate in a complicated multi-phase medium. 
The 24 /jm emission along the arc is likely to arise from 
the shocked molecular gas (e.g., OH and H 2 0), shock- 
heated dust grains, and even probably ions. The OH 
maser, HCO + , and H2O (if any) emissions are consistent 
with a C-type molecular shock. 

Second, let us discuss the scenario in which the arc 
represents the interstellar material swept up all the way 
by the SNR shock wave. The mass of the southeast- 
ern molecular arc is in the range of 0.5-1.1 x 10 4 e? 2 2 
Mq , with the gas mass observed in the whole line profile 
adopted as the upper limit and that of the broadened 
part as the lower limit (Table [2|). In the swept-up case, 
the southeastern arc consists of the molecular gas that 
was originally of mean density n(H 2 ) ~ 30-70(^2, cm -3 , 
distributed in an approximate 1 / 8th volume of a sphere 
of radius r s ~ 8'. 7 ~ 13ds.2 pc (considering the one- 
sided line broadening). In this scenario, the velocity 
of the SNR shock is represented by the arc's slow ex- 
pansion (v m ~ lOkrns -1 ), which would imply that the 
SNR is in the radiative phase. The explosion energy of 
the SNR, E~ 7 x 10 50 (r s /13pc) 312 [n(H 2 )/50cm- 3 ] 112 
(i! m /10kms -1 ) erg (Chevalier 1974), would have 
seemed to be normal, despite the seldom large age t ~ 
0.31r s /v m - 4 x 10 5 (r s /13pc) (« m /10 kms -1 ) -1 yr for 
X-ray-bright SNRs. However, the X-rays are bright along 
the southern shell and two X-ray emission peaks are es- 
sentially coincident with the two radio peaks (Fig. 02 to- 
ward which the pointing molecular line observation was 
made). A radiative shock as slow as 10 kins" 1 could 
not be responsible for such X-rays along the SNR rim. 
On the other hand, to accelerate particles to relativis- 
tic energies for emitting radio synchrotron, the lower 
limit of the shock velocity is (Draine & McKee 1993): 
v s > 64(x. i /10- 5 ) 1 / 8 [n(H 2 )/50cm- 3 ] 1 / 8 (^/O.l)- 1 / 4 

2 The mean density of the X-ray-emitting gas is estimated as 
~ 0.14ci 5 }J 2 cm -3 (where 6(5.2 = d/5.2kpc is used for scaling) from 
the volume emission measure n e n rrn H , cm 3 of 

an elliptical region (with half-axes 11'. 2 X 6'. 7), for which a volume 
of an oblate spheroid (with half-axes 11'. 2 X 11'. 2 X 6'. 7) is assumed. 



(T/100K) 01 kms _1 , where Xj is the ionization frac- 
tion, (j) CT the efficiency of the particle acceleration of 
the shock, and T the preshock molecular gas temper- 
ature. It is much larger than the observed expansion 
velocity of the molecular shell, and thus new parti- 
cle acceleration would be difficult to take place at this 
shock. The existing relativistic electrons would rapidly 
stream freely away from the shock front for a very low 
ionization fraction (Blandford & Cowie 1982). More- 
over, because of the ambipolar diffusion, the magnetic 
field will separate from the neutral gas in a timescale 
(Blandford & Cowie 1982): ~ 10 3 (<5r/l pc) 2 (xi/lQ- 5 ) 
[n(H 2 )/50cm- 3 ] 2 (B/10- 6 G)- 2 (s/100)~ 2 yr (where B 
is the magnetic field strength in the unshocked gas, s the 
density compression ratio, and 5r the shell thickness), 
very likely to be smaller than the remnant age in this 
case (even for s ~ 10). Therefore, such a slow expansion 
of molecular cloud is difficult to account for the SNR's 
radio emission. In view of the above points, this scenario 
does not apply to the molecular shell of Kes 69. 

Third, the molecular shell may be the debris of the 
cooled, condensed material, which was swept up by 
the stellar wind of the supernova progenitor from the 
molecular gas [n(H2) ~ 50cm~ 3 ] and is now hit by 
the SNR shock. This possibility could be compatible 
with the first scenario and naturally explain why there 
is pre-existing, very dense material that is now hit by 
the SNR shock. A molecular wind-bubble shell has re- 
cently been discovered coincident with the ring nebula 
G79.29+0.46 surrounding a luminous blue variable star 
(Rizzo et al. 2008), which provides an instance for this 
scenario. If there was a stellar- wind bubble, it was cre- 
ated (Castor et al. 1975; Weaver et al. 1977) ~ 1.4 x 

10 6 (r/13pc) 5 / 3 L3 6 1/3 [n(H 2 )/50cm^ 3 ] 1 / 3 yr ago, where 
L36 is the mechanical luminosity of the stellar wind in 
units of 10 36 ergs -1 , and the present velocity of the shell 
is - 6(r/13pc)~ 2 / 3 L^ 3 [n(H 2 )/50cm -3 ] -1 / 3 km s _1 . 
This velocity is comparable to the blueshift in the broad- 
ened line profile of the molecular gas along the shell. 

There appears to be a blowout morphology outlined 
by the extension of the radio/CO shell to the northeast 
out of the circle (Figs. 2] and [5]), somewhat similar to the 
blowout morphology in SNR N132D (e.g., Dickel & Milne 
1995; Xiao & Chen 2008). Actually, SNR N132D, which 
is in the vicinity of an MC, has been suggested to be 
shaped by the shock impacting on the stellar wind-bubble 
shell (Hughes 1987; Chen et al. 2003). The northeastern 
compact masers at both 69 and 85kms _1 are projected 
in the blowout region. In this case, as a possibility, the 
69kms~ 1 maser could arise from a dense clump deviat- 
ing from the systemic velocity by a strong perturbation, 
although it cannot be excluded from being nonassociated 
with the same SNR. 

4. SUMMARY 

We have performed a millimeter observation in 
CO and HCO+ lines toward SNR Kes 69. From the 
northeastern compact 1720 MHz OH maser region, the 
12 CO and 13 CO emission's peaks around 65kms~ 1 and 
85kms _1 , which are consistent with the masers' LSR 
velocities are detected. In the southeast, a molecular 
( 12 CO J=l-0) arc is revealed at 77-86 km s _1 , well 
coincident with the partial SNR shell detected in the 
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1.4 GHz radio continuum and mid-IR observations. 
An 85kms _1 HCO + emission is found to arise from a 
radio peak on the shell. Both the molecular arc and the 
HCO + emission at ~ SSkms^ 1 seem to be consistent 
with the presence of the extended 1720 MHz OH emis- 
sion along the southeastern boundary of Kes 69. The 
morphology correspondence between the CO emission 
and other band emission of the Kes 69 shell provides 
strong evidence for the association of the SNR with 
the ~ 85kms _1 component of molecular gas. There is 
another section of molecular arc at 79-86 km s" 1 in the 
northwest. Both the molecular arcs, together with the 
faint northern radio features, seem to be distributed 
along a circle of radius 8'. 7. The multiwavelength 
emissions along the southeastern shell can be explained 
by the impact of the SNR shock on a dense, clumpy 
patch of molecular gas. This pre-existing gas is likely to 
be a part of the cooled, clumpy debris of the interstellar 
molecular gas swept up by the progenitor's stellar wind. 
The association between SNR Kes 69 and the MC at 
the systemic velocity ~ 85kms _1 enables us to place 
the SNR at a kinematic distance of 5.2 kpc. 



Note added in proof. In another observation in CO 
lines toward SNR Kes 75, we similarly discover a molec- 
ular shell, a part of which follows the bright partial SNR 
shell seen in X-rays, mid-IR, and radio continuum, and 
provide effective evidence for the association between Kes 
75 and the adjacent MC (Su et al. 2009). Molecular 
shells are probably common in a number of SNRs but 
have rarely been studied. 
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